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Two-loop renormalization group equations in gauge theories with multiple U(l) groups are pre- 
sented. Instead of normalizing the abelian gauge fields in canonical forms, we retain kinetic-mixing 
terms and treat the mixing coefficients as free parameters. The /3- and 7-functions are then obtained 
in a straightforward manner. 
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I. INTRODUCTION 

Renormalization group equations (RGEs) relate pa- 
rameters at different scales. It plays an important role in 
the search of physics beyond the standard model (SM) 
and provides a unique window for physics at extremely 
high energy. Comprehensive study of two-loop RGEs 
in a general gauge theory had been performed long ago 
[1]. The calculation was first performed for a simple Lie 
group. By using suitable substitution rules, these results 
can be extended to cases of semi-simple Lie groups, as 
well as cases with only one U(l) sector. Recently, these 
calculations were updated by carefully treating the trans- 
formation properties of the fermion fields [2] and applied 
to the special case of the SM[3]. Extension to cases with 
multiple U(l) groups is more involved, due to the possible 
presence of kinetic mixing terms between different U(l) 
gauge bosons. In this letter, we extend the calculation to 
these cases. 

Gauge theories with multiple U(l) groups can arise in 
many context. In most grand unified theories, one in gen- 
eral has more than one U(l) sector. In string theories, 
one obtains extra U(l) sectors copiously. And there are 
no general rules to make these extra U(l) gauge bosons 
heavy. Were they discovered, extra U(l) bosons would 
provide important signatures of the underlying theory. 
In some model buildings, extra gauge U(l) bosons were 
sought to solve the fi problem in supersymmetric theories 
[4]. It is arguably that extra U(l) bosons are the most 
likely new physics next to supersymmetry to be discov- 
ered. Over the years, models with extra U(l)'s have been 
extensively studied [5]. 

Unique to multiple U(l) gauge theories, there could be 
kinetic mixing terms in the form 



T^my.v rpn 



(i) 



which are both Lorentz and gauge invariant. One needs 
extra symmetries to keep them out, while in non-abelian 
gauge theories, similar terms are forbidden automatically 
by gauge invariance. Even if one diagonalizes the the 
abelian gauge fields such that ^ m «(w 7^ n) vanish at 
one scale, the mixing terms will in general come back 
at another scale due to RGE runnings [6]. Especially, 
these mixings affect the various (3- and 7-functions, which 



should be included in consistent analysis. Convention- 
ally, one normalizes the abelian gauge fields into canon- 
ical forms, such that the kinetic mixing terms vanish at 
the tree level. For example, the /3-functions of gauge 
couplings have been obtained up to two-loops via this 
procedure in [7]. However, mixings usually come back 
in loop-levels and the whole procedure becomes rather 
complicated. In this letter, we will provide an alterna- 
tive approach. We retain the kinetic-mixing terms and 
treat £ mn as free parameters. It proves to be extremely 
convenient in practice. The (3- and 7-functions of the 
theory, including the /3-functions of the £ mn 's, can be ob- 
tained in a straightforward manner by extending existing 
calculations. 

The paper is organized as follows, section 2 presents 
our formalism and the 7-functions of the scalar and 
fermion fields. Section 3 presents the /3-functions of 
the gauge couplings and those of the £ mn 's. Section 4 
presents the /3-functions of the Yukawa couplings and sec- 
tion 5 those of quartic scalar couplings. We conclude in 
section 6. 



II. FORMALISM 

We start with a general renormalizable field theory 
with a group G x FJ m U m (l), where G is a compact sim- 
ple Lie group. The theory contains real scalar fields </> 
and two-component fermion fields ijij. The Lagrangian 
can be written as 



where 
Co 



C = Co + (gauge fixing + ghost terms) (2) 

* t TpmiMV rpn _ ^_ rpAfiu rpA 



h.c. 
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_ 2 {X^jC^k^a + h.c) - -^\abcd(t>a4'b4>c4'd, (3) 

where the diagonal terms with £ mm = 1 correspond to 
the usual kinetic terms for abelian gauge bosons, while 
the off-diagonal ^ m n(m ^ n) represent kinetic mixings. 
The gauge field strengths are 
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and covariant derivatives of matter fields are 

D^ a = d^ a -i 9 e*v£4 )b -iY,9m{Q s m )abV™4> b , 

m 
m 

where 9 A and t A are the representation matrices of G on 
the scalar and fermion fields, respectively, and and 

are their U m (l) charges. Since the scalar fields are 
real, A and are pure imaginary and antisymmetric. 

are real and diagonal. The constraint on Yukawa 
coupling matrices imposed by U(l) gauge invariance is 



Y\Q s m ) ba + Y a Q? n + Q? n Y a = 0, 



(4) 



which can be used to simplify the gauge structures. 

Conventionally, one normalizes the U(l) gauge fields 
to assume canonical forms before renormalization. It is 
easy to show that this normalization does not affect the 
/3-functions of g, Y a and X a bcd and the 7- functions of V A , 
4> a and tpj [8, 9]. On the other hand, the /3-functions of 
g rn in the conventional scheme are related to the ones in 
our scheme by a ^-dependent transformation. If one nor- 
malizes the U(l) gauge fields consistently, both schemes 
give the same physics prediction [9]. 

In the Landau gauge, the field propagators of the sim- 
ple gauge group G are 



D A ?(k)=5 AB l-g,„ + 



(5) 



(6) 



and the propagators of abelian gauge fields are 

n mn (k) - f- 1 (-a 1 k ^ kv \ 1 
^fiv V") — Wn I Ufiv ' ^2 J ^2 ' 

where is the inverse of £ mn . Note that, there are 
propagations from a U m (l) gauge field to a U n (l) gauge 
field due to kinetic mixings, with magnitude proportional 



For late convenience, we define 



nA n Bi 



Cf{S)=e A c 6t S 2 {S)5 AB = Trr 
Cf(F) = t A c t A , S 2 (F)6 AB = Tr[t A t B ], 
C 2 (G)S AB =f ACD f BCD , 
n2S -n s O s O 43 -O s O s O s O s 

2F -G F O F O af -Q F Q F Q F Q F 

Xmn = gmgn£, m m 



(7) 
(8) 
(9) 
(10) 

(11) 
(12) 



where in Eq. (12), there is no summation over m and n. 

Except for the £ m n(m ^ n), we separate the (3- function 
of each coupling constant x into two parts, 



0x = Px{G)+0 x (U), 



(13) 



where f3 x (G) is the result when one turns off all U(l) in- 
teractions and (3 X (U) is the additional contribution when 



the U(l) interactions are turned on. To two-loops, we 
have 



fix 
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Similarly, we separate the 7-function of field i 
1 
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(15) 



To two-loops, 13(G) and 7(G) have been given in [1, 2]. 
The additional contributions from U(l) groups, (3 X (U) 
and ji(U) are to be calculated in this letter. These can be 
obtained from existing calculations by closely inspecting 
the Feynman diagrams in [1] , and making proper substi- 
tution of gauge couplings and abelian gauge propagators. 
The process is straightforward but tedious. We content 
with by simply reporting the end results here. 



A. Scalar wave function renormalization 

To one loop, the additional contributions to the 7- 
functions of scalar fields induced by U(l) groups are 



(16) 



here and hereafter repeated indices such as m and n are 
implicitly summed over all abelian gauge groups, unless 
specified otherwise. To two loops, the additional contri- 
butions are 

lT\u) = x mn {5nTT[Qlf n (Y+ a Y b + Y+ b Y a )} 

+ ^Xpq(Qrnp)a b [imQl S q ) +40kTt(Q^)] 

+ \xpq(Q^ npq U + 3. 9 2 [c 2 (S)Ql s m ] ab ){n) 

where k = |(1) for two(four)-component fermions. 

B. Fermion wave function renormalization 

To one loop, the additional contributions to the 7- 
functions of fermion fields induced by abelian groups are 
zero. To two loops, 



lf{U) = Xr nn {-\Y a <?* n Y +a 
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\xpq{QZ)^Ql S q )+^{Q 2 r !q)] 

- lxpqQ 4 m F npq ~3g 2 C2(F)Q^ n }. (18) 



III. /3-FUNCTIONS OF GAUGE COUPLINGS 
AND KINETIC MIXINGS 

A. Nonabelian gauge coupling 
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At one loop level, the additional contributions to the 
/3-function of gauge coupling g of group G are zero. To 
two loop, one has 

pf \U) - 2g 3 Xmn T^ n (S) + 4Kg 3 Xmn T^ n (F) (19) 

where T^ n (S) and T^ n (F) are defined as 

T° n (S)5 AB = Tv(0 A e B Q 2 J n ) 
T% n {F)5 AB = Tr(t A t B QZ) 

B. Abelian gauge couplings 

The /3-function of the gauge coupling g p of group U p (l) 
is related to the anomalous dimension 7p of the corre- 
sponding gauge field V? [Figure 1(a)]. The Ward identity 
ensures that (3 P = g P j P . To two loops, 
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+ 2g 2 C 2 (S)S%(S)+WC2(F)S p 2 (F)} 



where 



Yf(F) = Tr [Cf(F)Y a Y +a ] , 



(20) 



(21) 



and there is no summation over p. C 2 (R) is the Casimir 
operator of G on the representations R. Cf(i?) and 
<Sf (-R) are the Casimir operator and the Dykin index of 
U p (l) acting on fields R, respectively. 



C. Kinetic mixings 

The /3-functions of the kinetic mixing coefficients £ pg 
are given by 



Ppq = {"fp + lq)ipq + ll 



(22) 



where repeated indices don't imply summation and 7^ is 
the anomalous dimension of the operator F viiv 'F? [Figure 
1(b)], 



7 « = 

Ipq 



9p9q 



3Tr(Q^) + -«Tr(Q^) 



(4tt) 2 

- ©{- ^ YaY+a QlO 

+ 4g 2 (Tr [Q 2 p S q C 2 (S)] + 2«Tr [Qf q C 2 (F)]) 
+ 4 Xm „ (Tr(Q 4 p S qmn ) + 2«Tr(Q^ mJ ) } (23 ) 
where the indices are summed over m(n) but not over 



FIG. 1: (a)Wave function renormalization of abelian gauge 
fields; (b)renormalization of abelian gauge kinetic mixing. 



IV. YUKAWA COUPLING 

The /3-functions of the Yukawa couplings can be ex- 
pressed as 



(3 a = 7 a + 7 +F Y a + Y a 7 F + 7f b Y 6 



(24) 



where 7° are the anomalous dimensions of the opera- 
tors (j) a ipj(ipk, J F and 7 f b are the anomalous dimensions 
of the corresponding fermions and bosons, respectively. 
To one loop level, the additional contributions to the f3- 
function of Yukawa couplings are 

fc\U) = -3 Xmn (Ql F n Y a + Y a QZ) (25) 
To two loop level, one has 

P { ?\U) = Xmn {3{Q 2 J n ,Y b Y+ a Y b } + 5Y b {QZ l ,Y+ a }Y b 
7 



[Q 2 J n Y b Y +b Y a + Y a Y +b Y b Q^ n 



- -(Y b Q™Y +b Y a + Y a Y+ b Q 2 ? n Y b ) 

+ &{Qf n Y a Y+ b Q F n Y b + Y b Qf n Y+ b Y a Q F n ) 

+ 5nY b Ti [Q 2 Z l {Y +a Y b + Y +b Y a )] 

+ 6 [(Q 2 m S n ) bc Y b Y +a Y c 2(Q 2 ^ n ) ac Y b Y+ c Y b ] 

+ \{Qmn)bc{Y b Y +c Y a + Y a Y+ c Y b ) 

2 2F 3 4F 

_ 3 S {QmvP^{F), Y a } - -^ Xpq {Q mnpq , Y a } 

+ 6g 2 C?(S){QZ,Y»} + 6g 2 (Q 2 ^ n ) ab {C 2 (F),Y b } 

+ SXpq(QmnU{Q 2 pq,Y b } 

+ -^X Pg {QZ> Y a }[llTr(Q 2 n s q ) + 40Tr(Q^)] 
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9 {Q mn C2{S))ab + 7,Xpq{Q mn p q )ab 



Y b 



" l X Pg (QZU[^(Qn F q )+ Tr (Qn S q )]Y b } (26) 



V. SCALAR QUARTIC COUPLING 

The /3-functions of the scalar quartic couplings can be 
expressed as 



Pabcd = labed + ^ "f S {i)^abcd 



(27) 
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where "/abed are the anomalous dimensions of the opera- 
tors 4>a4>b4 > c4'd, 7 S (i) is the anomalous dimension of the 
scalar field i. To one loop level, one has 



+ ZXmnX Pq AZZ q (28) 
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and to two loop level, one has 
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(29) 



The group factors are defined as 

Aatr = E c 2 Q ' m "« A ^ 



AfSr = E^(^)c 2 Q ' m "(z)A QbC(i 

y Q2F 5 ab = Tv [Q 2 J* n {Y +a Y b + Y +b Y a )] 

A ab7d P9 = 1 E X abef{Qt'Qp}ef{Qn,Qq}cd 

perms 



A x a £d n - \ E ^bef{QLo A }e f {Qle A } 

perms 

A ~ab7d Pq = 4 E A «^/{<3 m ,<5p}ce{Qf,Qg}d/ 

perms 



" perms 



r>Y,mnpq _ 
abed a 



\ Y.{QLQ s P }ab^(Q F n Q F q YCY+ d 

perms 

+ Y c Q F Q F q r +d ) 
B Y a Cd - \ E {Q m ,^} Q6 Tr(gr^r c r+ d 

perms 



perms 

Y c Q F t A Y +d " 



ZCd m = \ E {QLQ S P }abTr(Q F Y^Q F Y+ d ) 

perms 

= \ E {Q m ,^WTr(Q^Wy+ d ) 

perms 



rjY,mn 1 
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Q,mn 
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g (Qmn) fg^abef ^cdeg 

perms 

g X] ^abefKdgh(Qm)eg(Qn) 
perms 

J2c?> mn (i)H abcd 
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E Tr[{Q 2 rt F „,F }y+ 6 F c y+ d ] 



^abed 1 — g X {Om' Qp ^{Qfj }cd 

perms 

Amnpqkl \ A s-^Q .mn / -\ <\pqkl 

^abed - 2-^ °2 V l ) A abcd 

i 

A Q ab 7d n = E c 2 Q,mn (*)^^ 
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A Q,mnpq _ \^ /-iQ,mn/-\ *pq 
^abed — Z^° 2 y l )^abcd 
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A^J 1 = E^^d 
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A s a cr = E^(^™ 

i 

where C 2 Q ' mn (fc) are the eigenvalues of Q^Q^, and 
i . . . _ . _ 



■bed — g E {^ A '^ B }ab{^' 4 ' 

perms 

ff a6c d - E iTr(r a r+ b r c y+ d ) 

perms 



VI. CONCLUSION 



perms 



A Q £d m " P9 = E^'^^W^ 



We have now presented the two-loop RGEs for all di- 
mensionless parameters in a general gauge theory with 
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multiple abelian gauge groups. We have retained all ki- 
netic mixing terms and did not canonically normalize 
the abelian gauge fields. This approach proved to be 
extremely convenient. Were all abelian gauge interac- 
tions turned off, the two- loop RGEs were given in [1, 2]. 
When the abelian gauge interactions are turned on, the 
additional contributions have been obtained from exist- 
ing calculations by a close inspection of the relevant Feyn- 
man diagrams and by using suitable substitution rules. 
Were G not simple, but rather a direct product of simple 
groups Gk with gauge coupling g k , the RGEs can again 
be obtained by using substitution rules as given in [1, 2]. 
For part of g 4 -related terms, the following extra substi- 
tution rules are needed, 

i perms 

^fiT-k^E^W E {Gt,Q s m U{of,Q s n U 



g i s 2 {R)A^ d ^\gis k 2 (R) E {etQlU{e£,Q s n } cd 

perms 

g i C 2 (G)A" l b " d -> ~gtC% (G) ^ i^k . Qm}ab{^k^ Qn}cd 



(30) 



By introducing a dummy field, /3-functions of the dimen- 
sional parameters can also be deduced from those of Y a 
and \ abcd [2, 9]. 
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